SARS-CoV-2 escapes direct NK cell killing through Nsp1-mediated downregulation of ligands for NKG2D

SummaryNatural killer (NK) cells are cytotoxic effector cells that target and lyse virally-infected cells; many viruses therefore encode mechanisms to escape such NK cell killing. Here, we interrogated the ability of SARS-CoV-2 to modulate NK cell recognition and lysis of infected cells. We found that NK cells exhibit poor cytotoxic responses against SARS-CoV-2-infected targets, preferentially killing uninfected bystander cells. We demonstrate that this escape is driven by downregulation of ligands for the activating receptor NKG2D ("NKG2D-L"). Indeed, early in viral infection, prior to NKG2D-L downregulation, NK cells are able to target and kill infected cells; however, this ability is lost as viral proteins are expressed. Finally, we found that SARS-CoV-2 non-structural protein 1 (Nsp1) mediates downregulation of NKG2D-L and that Nsp1 alone is sufficient to confer resistance to NK cell killing. Collectively, our work reveals that SARS-CoV-2 evades NK cell cytotoxicity and describes a mechanism by which this occurs. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=190 SRC="FIGDIR/small/496341v1_ufig1.gif" ALT="Figure 1"> View larger version (47K): org.highwire.dtl.DTLVardef@bcffeeorg.highwire.dtl.DTLVardef@469b0eorg.highwire.dtl.DTLVardef@16dd205org.highwire.dtl.DTLVardef@f78070_HPS_FORMAT_FIGEXP M_FIG C_FIG

Activated interstitial macrophages are a predominant target of viral takeover and focus of inflammation in COVID-19 initiation in human lung

Early stages of deadly respiratory diseases such as COVID-19 have been challenging to elucidate due to lack of an experimental system that recapitulates the cellular and structural complexity of the human lung, while allowing precise control over disease initiation and systematic interrogation of molecular events at cellular resolution. Here we show healthy human lung slices cultured ex vivo can be productively infected with SARS-CoV-2, and the cellular tropism of the virus and its distinct and dynamic effects on host cell gene expression can be determined by single cell RNA sequencing and reconstruction of "infection pseudotime" for individual lung cell types. This revealed the prominent SARS-CoV-2 target is a population of activated interstitial macrophages, which as infection proceeds accumulate thousands of viral RNA molecules per cell, comprising up to 60% of the cellular transcriptome and including canonical and novel subgenomic RNAs. During viral takeover, there is cell-autonomous induction of a specific host interferon program and seven chemokines (CCL2, 7, 8, 13, CXCL10) and cytokines (IL6, IL10), distinct from the response of alveolar macrophages in which neither viral takeover nor induction of a substantial inflammatory response occurs. Using a recombinant SARS-CoV-2 Spike-pseudotyped lentivirus, we show that entry into purified human lung macrophages depends on Spike but is not blocked by cytochalasin D or by an ACE2-competing monoclonal antibody, indicating a phagocytosis- and ACE2-independent route of entry. These results provide a molecular characterization of the initiation of COVID-19 in human lung tissue, identify activated interstitial macrophages as a prominent site of viral takeover and focus of inflammation, and suggest targeting of these macrophages and their signals as a new therapeutic modality for COVID-19 pneumonia and progression to ARDS. Our approach can be generalized to define the initiation program and evaluate therapeutics for any human lung infection at cellular resolution.

Detection of bacterial co-infections and prediction of fatal outcomes in COVID-19 patients presenting to the emergency department using a 29 mRNA host response classifier.

ObjectiveClinicians in the emergency department (ED) face challenges in concurrently assessing patients with suspected COVID-19 infection, detecting bacterial co-infection, and determining illness severity since current practices require separate workflows. Here we explore the accuracy of the IMX-BVN-3/IMX-SEV-3 29 mRNA host response classifiers in simultaneously detecting SARS-CoV-2 infection, bacterial co-infections, and predicting clinical severity of COVID-19. Methods161 patients with PCR-confirmed COVID-19 (52.2% female, median age 50.0 years, 51% hospitalized, 5.6% deaths) were enrolled at the Stanford Hospital ED. RNA was extracted (2.5 mL whole blood in PAXgene Blood RNA) and 29 host mRNAs in response to the infection were quantified using Nanostring nCounter. ResultsThe IMX-BVN-3 classifier identified SARS-CoV-2 infection in 151 patients with a sensitivity of 93.8%. Six of 10 patients undetected by the classifier had positive COVID tests more than 9 days prior to enrolment and the remaining oscillated between positive and negative results in subsequent tests. The classifier also predicted that 6 (3.7%) patients had a bacterial co-infection. Clinical adjudication confirmed that 5/6 (83.3%) of the patients had bacterial infections, i.e. Clostridioides difficile colitis (n=1), urinary tract infection (n=1), and clinically diagnosed bacterial infections (n=3) for a specificity of 99.4%. 2/101 (2.8%) patients in the IMX-SEV-3 Low and 7/60 (11.7%) in the Moderate severity classifications died within thirty days of enrollment. ConclusionsIMX-BVN-3/IMX-SEV-3 classifiers accurately identified patients with COVID-19, bacterial co-infections, and predicted patients risk of death. A point-of-care version of these classifiers, under development, could improve ED patient management including more accurate treatment decisions and optimized resource utilization.

Comparative analysis of cell-cell communication at single-cell resolution

Inference of cell-cell communication (CCC) from single-cell RNA-sequencing data is a powerful technique to uncover putative axes of multicellular coordination, yet existing methods perform this analysis at the level of the cell type or cluster, discarding single-cell level information. Here we present Scriabin - a flexible and scalable framework for comparative analysis of CCC at single-cell resolution. We leverage multiple published datasets to show that Scriabin recovers expected CCC edges and use spatial transcriptomic data, genetic perturbation screens, and direct experimental manipulation of receptor-ligand interactions to validate that the recovered edges are biologically meaningful. We then apply Scriabin to uncover co-expressed programs of CCC from atlas-scale datasets, validating known communication pathways required for maintaining the intestinal stem cell niche and revealing species-specific communication pathways. Finally, we utilize single-cell communication networks calculated using Scriabin to follow communication pathways that operate between timepoints in longitudinal datasets, highlighting bystander cells as important initiators of inflammatory reactions in acute SARS-CoV-2 infection. Our approach represents a broadly applicable strategy to leverage single-cell resolution data maximally toward uncovering CCC circuitry and rich niche-phenotype relationships in health and disease.

TNFα-producing CD4+ T cells dominate the SARS-CoV-2-specific T cell response in COVID-19 outpatients and are associated with durable antibodies

SARS-CoV-2-specific CD4+ T cells are likely important in immunity against COVID-19, but our understanding of CD4+ longitudinal dynamics following infection and specific features that correlate with the maintenance of neutralizing antibodies remains limited. We characterized SARS-CoV-2-specific CD4+ T cells in a longitudinal cohort of 109 COVID-19 outpatients. The quality of the SARS-CoV-2-specific CD4+ response shifted from cells producing IFN{gamma} to TNF+ from five days to four months post-enrollment, with IFN{gamma}-IL21-TNF+ CD4+ T cells the predominant population detected at later timepoints. Greater percentages of IFN{gamma}-IL21-TNF+ CD4+ T cells on day 28 correlated with SARS-CoV-2 neutralizing antibodies measured seven months post-infection ({rho}=0.4, P=0.01). mRNA vaccination following SARS-CoV-2 infection boosted both IFN{gamma} and TNF producing, spike protein-specific CD4+ T cells. These data suggest that SARS-CoV-2-specific, TNF-producing CD4+ T cells may play an important role in antibody maintenance following COVID-19.

Magnitude and breadth of neutralizing antibody responses elicited by SARS-CoV-2 infection or vaccination

Multiple SARS-CoV-2 variants that possess mutations associated with increased transmission and antibody escape have arisen over the course of the current pandemic. While the current vaccines have largely been effective against past variants, the number of mutations found on the Omicron (B.1.529) spike appear to diminish the efficacy of pre-existing immunity. Using pseudoparticles expressing the spike of several SARS-CoV-2 variants, we evaluated the magnitude and breadth of the neutralizing antibody response over time in naturally infected and in mRNA-vaccinated individuals. We observed that while boosting increases the magnitude of the antibody response to wildtype (D614), Beta, Delta and Omicron variants, the Omicron variant was the most resistant to neutralization. We further observed that vaccinated healthy adults had robust and broad antibody responses while responses were relatively reduced in vaccinated pregnant women, underscoring the importance of learning how to maximize mRNA vaccine responses in pregnant populations. Findings from this study show substantial heterogeneity in the magnitude and breadth of responses after infection and mRNA vaccination and may support the addition of more conserved viral antigens to existing SARS-CoV-2 vaccines. One Sentence SummaryDiminished efficacy of pre-existing immunity to highly mutated SARS-CoV-2 variants.

SARS-CoV-2 infects human adipose tissue and elicits an inflammatory response consistent with severe COVID-19

The COVID-19 pandemic, caused by the viral pathogen SARS-CoV-2, has taken the lives of millions of individuals around the world. Obesity is associated with adverse COVID-19 outcomes, but the underlying mechanism is unknown. In this report, we demonstrate that human adipose tissue from multiple depots is permissive to SARS-CoV-2 infection and that infection elicits an inflammatory response, including the secretion of known inflammatory mediators of severe COVID-19. We identify two cellular targets of SARS-CoV-2 infection in adipose tissue: mature adipocytes and adipose tissue macrophages. Adipose tissue macrophage infection is largely restricted to a highly inflammatory subpopulation of macrophages, present at baseline, that is further activated in response to SARS-CoV-2 infection. Preadipocytes, while not infected, adopt a proinflammatory phenotype. We further demonstrate that SARS-CoV-2 RNA is detectable in adipocytes in COVID-19 autopsy cases and is associated with an inflammatory infiltrate. Collectively, our findings indicate that adipose tissue supports SARS-CoV-2 infection and pathogenic inflammation and may explain the link between obesity and severe COVID-19. One sentence summaryOur work provides the first in vivo evidence of SARS-CoV-2 infection in human adipose tissue and describes the associated inflammation.

Association Between SARS-CoV-2 RNAemia and Post-Acute Sequelae of COVID-19.

Determinants of Post-Acute Sequelae of COVID-19 are not known. Here we show that 75% of patients with viral RNA in blood (RNAemia) at presentation were symptomatic in the post-acute phase. RNAemia at presentation successfully predicted PASC, independent of patient demographics, initial disease severity, and length of symptoms.

Intramuscular SARS-CoV-2 vaccines elicit varying degrees of plasma and salivary antibody responses as compared to natural infection

Vaccination induced antibody and T-cell immune responses are important for systemic protection from COVID-19. Because SARS-CoV-2 infects and is transmitted by oral-pharyngeal mucosa, we wished to test mucosal antibodies elicited by natural infection or intramuscular vaccine injection. In a non-randomized observational study, we measured antibodies against the SARS-CoV-2 RBD in plasma and saliva from convalescent or vaccinated individuals and tested their neutralizing potential using a replication competent rVSV-eGFP-SARS-CoV-2. We found IgG and IgA anti-RBD antibodies as well as neutralizing activity in convalescent plasma and saliva. Two doses of mRNA vaccination (BNT162b2 or mRNA-1273) induced high levels of IgG anti-RBD in saliva, a subset of whom also had IgA, and significant neutralizing activity. We detected anti-RBD IgG and IgA with significant neutralizing potential in the plasma of single dose Ad26.COV2.S vaccinated individuals, and we detected slight amounts of anti-RBD antibodies in matched saliva. The role of salivary antibodies in protection against SARS-CoV-2 infection is unknown and merits further investigation. This study was not designed to, nor did it study the full kinetics of the antibody response or protection from infection, nor did it address variants of SARS-CoV-2.

Early immune responses have long-term associations with clinical, virologic, and immunologic outcomes in patients with COVID-19

The great majority of SARS-CoV-2 infections are mild and uncomplicated, but some individuals with initially mild COVID-19 progressively develop more severe symptoms. Furthermore, mild to moderate infections are an important contributor to ongoing transmission. There remains a critical need to identify host immune biomarkers predictive of clinical and virologic outcomes in SARS-CoV-2-infected patients. Leveraging longitudinal samples and data from a clinical trial of Peginterferon Lambda for treatment of SARS-CoV-2 infected outpatients, we used host proteomics and transcriptomics to characterize the trajectory of the immune response in COVID-19 patients within the first 2 weeks of symptom onset. We define early immune signatures, including plasma levels of RIG-I and the CCR2 ligands (MCP1, MCP2 and MCP3), associated with control of oropharyngeal viral load, the degree of symptom severity, and immune memory (including SARS-CoV-2-specific T cell responses and spike (S) protein-binding IgG levels). We found that individuals receiving BNT162b2 (Pfizer-BioNTech) vaccine had similar early immune trajectories to those observed in this natural infection cohort, including the induction of both inflammatory cytokines (e.g. MCP1) and negative immune regulators (e.g. TWEAK). Finally, we demonstrate that machine learning models using 8-10 plasma protein markers measured early within the course of infection are able to accurately predict symptom severity, T cell memory, and the antibody response post-infection.

The B.1.427/1.429 (epsilon) SARS-CoV-2 variants are more virulent than ancestral B.1 (614G) in Syrian hamsters

As novel SARS-CoV-2 variants continue to emerge, it is critical that their potential to cause severe disease and evade vaccine-induced immunity is rapidly assessed in humans and studied in animal models. In early January 2021, a novel variant of concern (VOC) designated B.1.429 comprising 2 lineages, B.1.427 and B.1.429, was originally detected in California (CA) and shown to enhance infectivity in vitro and decrease antibody neutralization by plasma from convalescent patients and vaccine recipients. Here we examine the virulence, transmissibility, and susceptibility to pre-existing immunity for B 1.427 and B 1.429 in the Syrian hamster model. We find that both strains exhibit enhanced virulence as measured by increased body weight loss compared to hamsters infected with ancestral B.1 (614G), with B.1.429 causing the most body weight loss among all 3 lineages. Faster dissemination from airways to parenchyma and more severe lung pathology at both early and late stages were also observed with B.1.429 infections relative to B.1. (614G) and B.1.427 infections. In addition, subgenomic viral RNA (sgRNA) levels were highest in oral swabs of hamsters infected with B.1.429, however sgRNA levels in lungs were similar in all three strains. This demonstrates that B.1.429 replicates to higher levels than ancestral B.1 (614G) or B.1.427 in the upper respiratory tract (URT) but not in the lungs. In multi-virus in-vivo competition experiments, we found that epsilon (B.1.427/B.1.429) and gamma (P.1) dramatically outcompete alpha (B.1.1.7), beta (B.1.351) and zeta (P.2) in the lungs. In the URT gamma, and epsilon dominate, but the highly infectious alpha variant also maintains a moderate size niche. We did not observe significant differences in airborne transmission efficiency among the B.1.427, B.1.429 and ancestral B.1 (614G) variants in hamsters. These results demonstrate enhanced virulence and high relative fitness of the epsilon (B.1.427/B.1.429) variant in Syrian hamsters compared to an ancestral B.1 (614G) strain. Author SummaryIn the last 12 months new variants of SARS-CoV-2 have arisen in the UK, South Africa, Brazil, India, and California. New SARS-CoV-2 variants will continue to emerge for the foreseeable future in the human population and the potential for these new variants to produce severe disease and evade vaccines needs to be understood. In this study, we used the hamster model to determine the epsilon (B.1.427/429) SARS-CoV-2 strains that emerged in California in late 2020 cause more severe disease and infected hamsters have higher viral loads in the upper respiratory tract compared to the prior B.1 (614G) strain. These findings are consistent with human clinical data and help explain the emergence and rapid spread of this strain in early 2021.

An intranasal ASO therapeutic targeting SARS-CoV-2

The COVID-19 pandemic is exacting an increasing toll worldwide, with new SARS-CoV-2 variants emerging that exhibit higher infectivity rates and that may partially evade vaccine and antibody immunity1. Rapid deployment of non-invasive therapeutic avenues capable of preventing infection by all SARS-CoV-2 variants could complement current vaccination efforts and help turn the tide on the COVID-19 pandemic2. Here, we describe a novel therapeutic strategy targeting the SARS-CoV-2 RNA using locked nucleic acid antisense oligonucleotides (LNA ASOs). We identified an LNA ASO binding to the 5 leader sequence of SARS-CoV-2 ORF1a/b that disrupts a highly conserved stem-loop structure with nanomolar efficacy in preventing viral replication in human cells. Daily intranasal administration of this LNA ASO in the K18-hACE2 humanized COVID-19 mouse model potently (98-99%) suppressed viral replication in the lungs of infected mice, revealing strong prophylactic and treatment effects. We found that the LNA ASO also represses viral infection in golden Syrian hamsters, and is highly efficacious in countering all SARS-CoV-2 "variants of concern" tested in vitro and in vivo, including B.1.427, B.1.1.7, and B.1.351 variants3. Hence, inhaled LNA ASOs targeting SARS-CoV-2 represents a promising therapeutic approach to reduce transmission of variants partially resistant to vaccines and monoclonal antibodies, and could be deployed intranasally for prophylaxis or via lung delivery by nebulizer to decrease severity of COVID-19 in infected individuals. LNA ASOs are chemically stable and can be flexibly modified to target different viral RNA sequences4, and they may have particular impact in areas where vaccine distribution is a challenge, and could be stockpiled for future coronavirus pandemics.

SARS-CoV-2 subgenomic RNA kinetics in longitudinal clinical samples

BackgroundGiven the persistence of viral RNA in clinically recovered COVID-19 patients, subgenomic RNAs (sgRNA) have been reported as potential molecular viability markers for SARS-CoV-2. However, few data are available on their longitudinal kinetics, compared with genomic RNA (gRNA), in clinical samples. MethodsWe analyzed 536 samples from 205 patients with COVID-19 from placebo-controlled, outpatient trials of Peginterferon Lambda-1a (Lambda; n=177) and favipiravir (n=359). Nasal swabs were collected at three time points in the Lambda (Day 1, 4 and 6) and favipiravir (Day 1, 5, and 10) trials. N-gene gRNA and sgRNA were quantified by RT-qPCR. To investigate the decay kinetics in vitro, we measured gRNA and sgRNA in A549ACE2+ cells infected with SARS-CoV-2, following treatment with remdesivir or DMSO control. ResultsAt six days in the Lambda trial and ten days in the favipiravir trial, sgRNA remained detectable in 51.6% (32/62) and 49.5% (51/106) of the samples, respectively. Cycle threshold (Ct) values for gRNA and sgRNA were highly linearly correlated (Pearsons r=0.87) and the rate of increase did not differ significantly in Lambda (1.36 cycles/day vs 1.36 cycles/day; p = 0.97) or favipiravir (1.03 cycles/day vs 0.94 cycles/day; p=0.26) trials. From samples collected 15-21 days after symptom onset, sgRNA was detectable in 48.1% (40/83) of participants. In SARS-CoV-2 infected A549ACE2+ cells treated with remdesivir, the rate of Ct increase did not differ between gRNA and sgRNA. ConclusionsIn clinical samples and in vitro, sgRNA was highly correlated with gRNA and did not demonstrate different decay patterns to support its application as a viability marker. SummaryWe observed prolonged detection of subgenomic RNA in nasal swabs and equivalent decay rates to genomic RNA in both longitudinal nasal swabs and in remdesivir-treated A549ACE2+ cells infected with SARS-CoV-2. Taken together, these findings suggest that subgenomic RNA from SARS-CoV-2 is comparably stable to genomic RNA and that its detection is therefore not a more reliable indicator of replicating virus.

Multi-omic profiling reveals widespread dysregulation of innate immunity and hematopoiesis in COVID-19

Our understanding of protective vs. pathologic immune responses to SARS-CoV-2, the virus that causes Coronavirus disease 2019 (COVID-19), is limited by inadequate profiling of patients at the extremes of the disease severity spectrum. Here, we performed multi-omic single-cell immune profiling of 64 COVID-19 patients across the full range of disease severity, from outpatients with mild disease to fatal cases. Our transcriptomic, epigenomic, and proteomic analyses reveal widespread dysfunction of peripheral innate immunity in severe and fatal COVID-19, with the most profound disturbances including a prominent neutrophil hyperactivation signature and monocytes with anti-inflammatory features. We further demonstrate that emergency myelopoiesis is a prominent feature of fatal COVID-19. Collectively, our results reveal disease severity-associated immune phenotypes in COVID-19 and identify pathogenesis-associated pathways that are potential targets for therapeutic intervention. One Sentence SummarySingle-cell profiling demonstrates multifarious dysregulation of innate immune phenotype associated with COVID-19 severity.

SARS-CoV-2 RNAaemia predicts clinical deterioration and extrapulmonary complications from COVID-19.

BackgroundThe determinants of COVID-19 disease severity and extrapulmonary complications (EPCs) are poorly understood. We characterise the relationships between SARS-CoV-2 RNAaemia and disease severity, clinical deterioration, and specific EPCs. MethodsWe used quantitative (qPCR) and digital (dPCR) PCR to quantify SARS-CoV-2 RNA from nasopharyngeal swabs and plasma in 191 patients presenting to the Emergency Department (ED) with COVID-19. We recorded patient symptoms, laboratory markers, and clinical outcomes, with a focus on oxygen requirements over time. We collected longitudinal plasma samples from a subset of patients. We characterised the role of RNAaemia in predicting clinical severity and EPCs using elastic net regression. Findings23{middle dot}0% (44/191) of SARS-CoV-2 positive patients had viral RNA detected in plasma by dPCR, compared to 1{middle dot}4% (2/147) by qPCR. Most patients with serial measurements had undetectable RNAaemia 10 days after onset of symptoms, but took 16 days to reach maximum severity, and 33 days for symptoms to resolve. Initially RNAaemic patients were more likely to manifest severe disease (OR 6{middle dot}72 [95% CI, 2{middle dot}45 - 19{middle dot}79]), worsening of disease severity (OR 2{middle dot}43 [95% CI, 1{middle dot}07 - 5{middle dot}38]), and EPCs (OR 2{middle dot}81 [95% CI, 1{middle dot}26 - 6{middle dot}36]). RNA load correlated with maximum severity (r = 0{middle dot}47 [95% CI, 0{middle dot}20 - 0{middle dot}67]). InterpretationdPCR is more sensitive than qPCR for the detection of SARS-CoV-2 RNAaemia, which is a robust predictor of eventual COVID-19 severity and oxygen requirements, as well as EPCs. Since many COVID-19 therapies are initiated on the basis of oxygen requirements, RNAaemia on presentation might serve to direct early initiation of appropriate therapies for the patients most likely to deteriorate. FundingNIH/NIAID (Grants R01A153133, R01AI137272, and 3U19AI057229 - 17W1 COVID SUPP #2) and a donation from Eva Grove. Research in contextO_ST_ABSEvidence before this studyC_ST_ABSThe varied clinical manifestations of COVID-19 have directed attention to the distribution of SARS-CoV-2 in the body. Although most concentrated and tested for in the nasopharynx, SARS-CoV-2 RNA has been found in blood, stool, and numerous tissues, raising questions about dissemination of viral RNA throughout the body, and the role of this process in disease severity and extrapulmonary complications. Recent studies have detected low levels of SARS-CoV-2 RNA in blood using either quantitative reverse transcriptase real-time PCR (qPCR) or droplet digital PCR (dPCR), and have associated RNAaemia with disease severity and biomarkers of dysregulated immune response. Added value of this studyWe quantified SARS-CoV-2 RNA in the nasopharynx and plasma of patients presenting to the Emergency Department with COVID-19, and found an array-based dPCR platform to be markedly more sensitive than qPCR for detection of SARS-CoV-2 RNA, with a simplified workflow well-suited to clinical adoption. We collected serial plasma samples during patients course of illness, and showed that SARS-CoV-2 RNAaemia peaks early, while clinical condition often continues to worsen. Our findings confirm the association between RNAaemia and disease severity, and additionally demonstrate a role for RNAaemia in predicting future deterioration and specific extrapulmonary complications. Implications of all the available evidenceVariation in SARS-CoV-2 RNAaemia may help explain disparities in disease severity and extrapulmonary complications from COVID-19. Testing for RNAaemia with dPCR early in the course of illness may help guide patient triage and management.

Additional analyses exploring the hypothesized transdifferentiation of plasmablasts to developing neutrophils in severe COVID-19

We thank Alquicira-Hernandez et al. for their reanalysis of our single-cell transcriptomic dataset profiling peripheral immune responses to severe COVID-19. We agree that careful analysis of single-cell sequencing data is important for generating cogent hypotheses but find several aspects of their criticism of our analysis to be problematic. Here we respond briefly to misunderstandings and inaccuracies in their commentary that may have led to misinformed interpretation of our results.

Integrated analysis of multimodal single-cell data

The simultaneous measurement of multiple modalities, known as multimodal analysis, represents an exciting frontier for single-cell genomics and necessitates new computational methods that can define cellular states based on multiple data types. Here, we introduce weighted-nearest neighbor analysis, an unsupervised framework to learn the relative utility of each data type in each cell, enabling an integrative analysis of multiple modalities. We apply our procedure to a CITE-seq dataset of hundreds of thousands of human white blood cells alongside a panel of 228 antibodies to construct a multimodal reference atlas of the circulating immune system. We demonstrate that integrative analysis substantially improves our ability to resolve cell states and validate the presence of previously unreported lymphoid subpopulations. Moreover, we demonstrate how to leverage this reference to rapidly map new datasets, and to interpret immune responses to vaccination and COVID-19. Our approach represents a broadly applicable strategy to analyze single-cell multimodal datasets, including paired measurements of RNA and chromatin state, and to look beyond the transcriptome towards a unified and multimodal definition of cellular identity. AvailabilityInstallation instructions, documentation, tutorials, and CITE-seq datasets are available at http://www.satijalab.org/seurat

SARS-CoV-2 Antibody Responses Correlate with Resolution of RNAemia But Are Short-Lived in Patients with Mild Illness

SARS-CoV-2-specific antibodies, particularly those preventing viral spike receptor binding domain (RBD) interaction with host angiotensin-converting enzyme 2 (ACE2) receptor, could offer protective immunity, and may affect clinical outcomes of COVID-19 patients. We analyzed 625 serial plasma samples from 40 hospitalized COVID-19 patients and 170 SARS-CoV-2-infected outpatients and asymptomatic individuals. Severely ill patients developed significantly higher SARS-CoV-2-specific antibody responses than outpatients and asymptomatic individuals. The development of plasma antibodies was correlated with decreases in viral RNAemia, consistent with potential humoral immune clearance of virus. Using a novel competition ELISA, we detected antibodies blocking RBD-ACE2 interactions in 68% of inpatients and 40% of outpatients tested. Cross-reactive antibodies recognizing SARS-CoV RBD were found almost exclusively in hospitalized patients. Outpatient and asymptomatic individuals serological responses to SARS-CoV-2 decreased within 2 months, suggesting that humoral protection may be short-lived.

Progenitor identification and SARS-CoV-2 infection in long-term human distal lung organoid cultures

The distal lung contains terminal bronchioles and alveoli that facilitate gas exchange and is affected by disorders including interstitial lung disease, cancer, and SARS-CoV-2-associated COVID-19 pneumonia. Investigations of these localized pathologies have been hindered by a lack of 3D in vitro human distal lung culture systems. Further, human distal lung stem cell identification has been impaired by quiescence, anatomic divergence from mouse and lack of lineage tracing and clonogenic culture. Here, we developed robust feeder-free, chemically-defined culture of distal human lung progenitors as organoids derived clonally from single adult human alveolar epithelial type II (AT2) or KRT5+ basal cells. AT2 organoids exhibited AT1 transdifferentiation potential, while basal cell organoids progressively developed lumens lined by differentiated club and ciliated cells. Organoids consisting solely of club cells were not observed. Upon single cell RNA-sequencing (scRNA-seq), alveolar organoids were composed of proliferative AT2 cells; however, basal organoid KRT5+ cells contained a distinct ITGA6+ITGB4+ mitotic population whose proliferation segregated to a TNFRSF12Ahi subfraction. Clonogenic organoid growth was markedly enriched within the TNFRSF12Ahi subset of FACS-purified ITGA6+ITGB4+ basal cells from human lung or derivative organoids. In vivo, TNFRSF12A+ cells comprised ~10% of KRT5+ basal cells and resided in clusters within terminal bronchioles. To model COVID-19 distal lung disease, we everted the polarity of basal and alveolar organoids to rapidly relocate differentiated club and ciliated cells from the organoid lumen to the exterior surface, thus displaying the SARS-CoV-2 receptor ACE2 on the outwardly-facing apical aspect. Accordingly, basal and AT2 "apical-out" organoids were infected by SARS-CoV-2, identifying club cells as a novel target population. This long-term, feeder-free organoid culture of human distal lung alveolar and basal stem cells, coupled with single cell analysis, identifies unsuspected basal cell functional heterogeneity and exemplifies progenitor identification within a slowly proliferating human tissue. Further, our studies establish a facile in vitro organoid model for human distal lung infectious diseases including COVID-19-associated pneumonia.

Human B cell clonal expansion and convergent antibody responses to SARS-CoV-2

During virus infection B cells are critical for the production of antibodies and protective immunity. Here we show that the human B cell compartment in patients with diagnostically confirmed SARS-CoV-2 and clinical COVID-19 is rapidly altered with the early recruitment of B cells expressing a limited subset of IGHV genes, progressing to a highly polyclonal response of B cells with broader IGHV gene usage and extensive class switching to IgG and IgA subclasses with limited somatic hypermutation in the initial weeks of infection. We identify extensive convergence of antibody sequences across SARS-CoV-2 patients, highlighting stereotyped naive responses to this virus. Notably, sequence-based detection in COVID-19 patients of convergent B cell clonotypes previously reported in SARS-CoV infection predicts the presence of SARS-CoV/SARS-CoV-2 cross-reactive antibody titers specific for the receptor-binding domain. These findings offer molecular insights into shared features of human B cell responses to SARS-CoV-2 and other zoonotic spillover coronaviruses.